Dephasing dynamics of Rydberg atom spin waves 
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A theory of Rydberg atom interactions is used to derive analytical forms for the spin wave pair 
correlation function in laser-excited cold-atom vapors. This function controls the quantum statistics 
of light emission from dense, inhomogeneous clouds of cold atoms of various spatial dimensionalities. 
The results yield distinctive scaling behaviors on the microsecond timescale, including generalized 
exponential decay. A detailed comparison is presented with a recent experiment on a cigar-shaped 
atomic ensemble [Y. Dudin and A. Kuzmich, Science 336, 887 (2012)], in which Rb atoms are 
excited to a set of Rydberg levels. 



Dipolc-dipole interactions between atoms excited to 
Rydberg levels are very strong compared to interactions 
between atoms in their ground states A novel ap- 
plication, proposed in Ref. [2|, is that the Rydberg in- 
teraction could be utilized to entangle atoms, even sep- 
arated by up to 10 microns, by controlling the atomic 
phase shifts accumulated in the presence of laser excita- 
tion. As a consequence Rydberg interactions offer great 
promise — e.g., in the areas of ultracold atomic physics 
and quantum information physics — for the generation 
and control of many-body quantum states [3J, sources 
of nonclassical light emission and the realization of fast 
quantum gates | ]j . 

The Rydberg blockade is a paradigm for strong Ryd- 
berg atom interactions: an atom, once excited to a Ry- 
dberg level, induces a level shift of nearby atoms that 
prevents their excitation Q . This process promises near- 
perfect efficiency for single photon emission and realiza- 
tion of quantum gates. For the blockade to be effective, 
the interaction strength between every atom pair should 
exceed the Rabi frequency and linewidth of the excita- 
tion laser; this requirement imposes stringent limits on 
the size of the atomic cloud. 

A different interaction regime may be reached by load- 
ing an atomic ensemble with a Poissonian number dis- 
tribution of laser-induced Rydberg excitations [6j. In- 
teractions then cause temporal dephasing of states with 
more than one excitation: an effect based on the varia- 
tion of the interatomic potential in a spatially extended 
ensemble. The retrieval process maps the spin wave pair 
correlations onto the emitted light, via a phase match- 
ing condition |7J, promising a fast, high-quality, single- 
photon source. The efficiency of the source is limited by 
the probability amplitude of the single-excitation com- 
ponent of the many-body wavefunction. The experimen- 
tal observation of spin wave dephasing has recently been 
reported [8j]. That the interaction mechanism is dom- 
inated by dephasing, as opposed to Rydberg blockade, 
was inferred from the evolution of the excited fraction of 
Rydberg atoms as a function of the Rabi frequency and 
by the absence of many-body oscillations, typical of the 
blockade regime. 

In this Paper, motivated by the desire to obtain a 
closed expression for the time-dependence of the spin 
wave pair correlation function, we present a theory of the 



spin wave dynamics based on Rydberg atom short- and 
long-ranged interactions. By using asymptotic methods 
appropriate for the long-time regime, we obtain analyt- 
ical expressions for the dynamics of the spin wave pair 
correlation function, which controls the sub-Poissonian 
quantum statistics of the emitted light field and the speed 
of the single-photon protocol. In order to be able to ad- 
dress highly asymmetric clouds, we derive expressions for 
samples of varying spatial dimensionality D. We present 
results for both van der Waals (vdW) and dipole-dipole 
interaction potentials, and show that experimental data 
reported in Ref. [|[ is in agreement with the vdW theory. 

Spin wave dephasing. We consider a gas of iV atoms of 
fixed positions r M (with /i = 1, . . . , N), laser excited from 
the atomic ground state \g) to a Rydberg level \r). Wc 
denote by \G) the JV-atom ground state, and we denote 
by \fiis) = \g±, . . . , r^, . . . , rv, . . . ,gjv) the state in which 
atoms y, and v (and only those atoms) are excited from 
\g) to |r). In the limit in which only a few atoms are ex- 
cited, we can effectively describe collective excitations by 
means of quasi-bosonic spin waves associated with the de- 
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where ko is the wavevector of the spin wave, which is 
fixed by the excitation process. We assume that the laser 
excitation prepares the system at time t = in the state 
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As the excitation is assumed to be weak, the distribution 
|c m | 2 is typically a non- increasing function of m. Two- 
body interactions between Rydberg-excited atoms pL and 
v give a level-shift A M!/ = V^/H that depends on the spe- 
cific form of the interaction potential V^ v . Here, we shall 
assume that V^„ = C Q /|r M — rv| Q , with a = 3 or 6 and 
C a determined by the Rydberg target state 0, 0] . After 
excitation, the ./V atoms interact for a time t. As a result 
the atomic states are phase shifted; e.g., in the two-body 
sector, \fj,u) — > exp(— iA^t^/iv) . More generally, the 
phase shift for a state with m excitations, \fii . . . fj, m ) , 
is given by ^...^ = -i£i<i<i<m A M<W . The spin- 
wave pair correlation function, after storage time t, 
gW{t) ee (^ W^ (t)^ ko (t)4 o (i))/(^ o (t)^ ko (t))2, is 
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*iven, with j m -k = m!/[(m — k)\N m+k ], by 
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where we have assumed that A^ 3> 1 . The term for m = 1 
in the denominator corresponds to |ci| 2 , as states having 
a single excitation are not phase shifted. For m = 2 in 
the numerator, we have only the sum within the abso- 
lute value that we define as the atom pair interference 
function 
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It is known that g^(0) = 1 for a Poisson distribution 
of c m . Interactions cause destructive interference of the 
distinct atom pair contributions, resulting in temporal 
decay of V(t). Higher-order interference contributions 
may be expressed in terms of the two-body result 
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where we use the results of Ref. [10J. For long enough 
times, t > t c where t c is a characteristic time associated 
with the specific interaction potential, \V(t)\ <C 1 and, 
therefore, the correlation function reduces to the form [6j 
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In Ref. we exhibit the direct relationship between 
states of spin wave excitation and those of photons emit- 
ted into a phase-matched mode of the field. As a conse- 
quence, the spin wave pair correlation function is equal to 
the Glauber's normalized second-order correlation func- 
tion of the phase-matched field mode. 
Dephasing in a cigar-shaped ensemble. In this section, 
by using formula ©, we numerically calculate g^ 2 \t) for 
a cigar-shaped ensemble and compare the results with 
the data of the experiment described in Ref. @. In that 
work, the spin wave correlations were inferred by mea- 
suring the probability of two coincident photoelectric de- 
tection events, in a phase matched mode determined by 
the experimental setup. Instead of studying the tem- 
poral dynamics for a given Rydberg level, the principal 
quantum number n was varied for a fixed storage-time 
T s = 0.3 /us. As the phase shifts are given by V^ v T s /h, 
and the Rydberg interaction scales with n, the measure- 
ments effectively map out the dynamics. The longitu- 
dinal size of the ensemble is set by the optical lattice 
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Figure 1. (Color online). Behavior of the spin wave pair 
correlation function (?' 2 ' vs. n for an asymmetrical cloud hav- 
ing (w z , w±) — (15, 6.4) /tm, compared with the experimental 
data from Ref. [g] (Black dots). Also shown are results from 
the double spin wave model (Blue, dashed line). Experimental 
data are reprinted with permission from AAAS and authors. 



beam, of waist w z = 15/im; the transverse area is de- 
termined by the overlap of two laser beams, each of 
waist 9/^m, giving a transverse waist of the excitation 
zone w± of (9/v2) /^m. For an ensemble of peak den- 
sity po ~ 10 12 cm -3 , the number of interacting atoms is 
given by N = (7r/2) 3 / 2 /9 w z w\ ~ 1200. In the present 
work, we determine the interatomic potentials for dif- 
ferent n via a single-channel model, \ns\/2, fiSi/2) ^ 
\ n P3/2i i n ~ 1)^3/2)) which is based on a semiclassical cal- 
culation of the dipole transition matrix elements Q . We 
evaluate the atom pair interference function by Monte 
Carlo integration, and use it to determine the contribu- 
tion of multiple excitations: the result converges once we 
have included of the order of 10 excitations. We sam- 
ple A" = 1000 atoms, although the result is essentially 
independent of N, provided A~ is much larger than the 
number of excitations. We assume that the initial spin- 
wave distribution state is Poissonian ([1} with unit mean, 
i.e., c m = lA/em!. To facilitate comparison with the 
data of Ref. the theoretical values are adjusted to ac- 
count for the measured photodetection background noise 



(1 — 9bq)9^ + 9bq ■ Figure [T] shows striking 
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agreement between the model and the experimental data 
over the full range of principal quantum numbers mea- 
sured. We note that a model retaining only single and 
double excitations is good for large principal quantum 
numbers, i.e., n > 60. 

Analytical treatment of spin wave dephasing. To develop 
an analytical theory of the spin wave dynamics in the 
limit AT 3> 1, we define V(t) via 
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where p(A) is the probability distribution for the 
interaction-induced frequency shift A associated with the 
ensemble of atom locations. We note that p(A) is deter- 
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mined by the interaction potential and by the geometry 
of the atomic ensemble. 

To determine p, we consider a 13-dimensional isotropic, 
ensemble having a Gaussian density profile; the formal- 
ism makes it possible to deal qualitatively also with 
anisotropic samples whose aspect ratio can be experi- 
mentally engineered via a combination of trap and excita- 
tion beams [3j . The two-particle spatial distribution P is 
given by Pd{yi,y 2 ) = Nd e~( r i +r 2)/( 2<T ), where w = 2a 
is the waist of the cloud, and the normalization Md is 
given by Nd = {2ira 2 )~ D . Then the probability den- 
sity for finding two particles separated by the vector R 
is given by 

Pd(R) = J d D n d D r 2 P{v 1 ,v 2 )5{H-Y 1 +v 2 ) 



= ( 47r(T 2)£>/2- 

Hence, one has that the 
particle separations R = 

by P D (R) = S D -ie- R 'K 
[= {2n D / 2 /T{D/2))R D - 1 
area of the surface of a D 
From this distribution, we 
distribution of a function 
bian transformation: 
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distribution Pd{R) of inter- 
|R| in dimension D is given 

4CT2 V(47ra 2 ) D / 2 , where S D ^ 
] is the (D — l)-dimensional 
-dimensional ball of radius R. 
have access to the probability 
f(R), via the standard Jaco- 
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We use this result to determine the probability distribu- 
tion of Rydberg atom level shifts for different interaction 
potentials, enabling calculation of the atom pair interfer- 
ence function via asymptotic methods. 
"Isotropic" dipole-dipole interaction. The interaction be- 
tween states of nonzero dipole moment is proportional to 
C3 / R 3 , where C3 is determined by the contributing scat- 
tering channels [9] . One way to achieve this interaction is 
by coupling near-degenerate atomic orbitals of different 
parity via an electrostatic field, so that the atoms acquire 
a dipole moment aligned along the field . For systems 
of reduced dimensionality (i.e., D being 1 or 2), the in- 
teraction can be made isotropic by applying the external 
field in a direction orthogonal to the system; for the case 
of D = 3, isotropy is not truly achievable, so an isotropic 
theory should provide only a qualitative description. 

We now calculate the atom pair interference function: 
V${t) = r o °° dupffi^e-™^ 1 . Here, we define the di- 
mensionless time r = flddt, where hfldd = (C^/a 3 ), and 
the normalized interaction energy k = (a/R) 3 , whose 
probability distribution is then given by Eq. ([8]l. 
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By using the method of steepest descents (see, e.g., 
Ref. [IH), we evaluate the long-time asymptotic behavior 
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Figure 2. (Color online). Asymptotic behavior of atom pair 
interference function for double excitations of the 100s level, 
(a) Isotropic dipole-dipole interactions: solid lines are numer- 
ical simulation for D = 1 [blue (black in BW)], 2 [red (dark 
gray)] and 3 [green (light gray)], compared with the long-time 
asymptotic behavior given by the modulus square of Eq. (|10|) 
(dashed lines). The green (light gray) dot-dashed line includes 
the first asymptotic correction for D — 3. (b) vdW interac- 
tions: (color legend as in Fig. a) the long-time behavior shows 
the square modulus of formula (|12l) . The first-order asymp- 
totic correction for D = 3 is indistinguishable from the solid 
line. Parameters are given in the text. 



of the integral, obtaining at leading order 
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where B = (5/2) x 6~ 3 / 5 e""/ 5 , and A D = y/n/5 x 
e -<*(D-i)/io 6 (D-i)/6/(2 D - 2 r(I>/2)). According to 
Eq. ([5]), the asymptotic decay of follows the square 



modulus of Vj)(t), which is given by 



\r d n d (r)\ 2 ^ 



D 

\A D \ 2 t I(d-i) c -2R C [b]t 2 '\ We note that thc next 
asymptotic correction to Vjj{t), multiplies the latter by 
the complex factor {1 + e'f /(6t) 2 / 5 [§(1 + f)(D - 2) + 
i|]}. In Fig. [2^, discussed further below, we compare 
the asymptotics with numerical results generated for an 
ensemble of atoms at randomly distributed positions [l3[ ■ 
Van der Waals interaction. In the absence of external 
electromagnetic fields, the interaction between atoms ex- 
cited to Rydberg levels has a dipole-dipole character at 
small separation and crosses over to the vdW (i.e., i?~ 6 ) 
form at larger separations As the typical size of 
the ensembles under consideration is much larger than 
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the cross-over radius, we assume a vdW interaction over 
the whole cloud. In analogy with the discussion of the 
dipole-dipole case, we define a dimensionless vdW energy 
yc = (er/i?) 6 , distributed according to 
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An asymptotic expression for the atom pair interference 
function © is found using the method of steepest de- 
scents as before, 
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dt 8 e" 



Ft 1 '* 



(12) 



where the dimensionless time is r = Q v dwt> 
Kl v dw = (Ce/c 6 ), and we define the coefficients 
C D = V2^e- i ra(£ , -i)i2( £> - 1 )/ 8 / [2 £> r(£>/2)] and T = 
(l2 1 / 4 /3) e _llr / 8 . We note that the next order asymp- 
totic correction maps V^ dw (t) V v D dw {t) x {1 + 

e-/8 /(12T )V4[3 (1 + D ){ SD _ 1)+ 36] }< In Fig . we 

show the dephasing dynamics for the vdW interaction for 
atomic ensembles of various dimensionality but common 
waist size. 

These results may be compared with those for dipole- 
dipole interactions, shown in Fig. [2^. We compare the 
asymptotics with numerical calculations for the target 
Rydberg level 100s, for which C 3 /H = 2tt x 10 5 MHz /im 3 
and Ca/h = 2n x 5.3 x 10 7 MHz^im 6 . For an ensem- 
ble of waist w = 30 /im, we obtain typical level shifts 
of n dd = 2tt x 30 MHz and fl vdw = 2n x 4.7 MHz. We 
note that for the vdW case, decay is roughly two or- 
ders of magnitude slower than for the dipole-dipole in- 
teraction. This result is dependent on the size of the 
ensemble, which sets the typical interaction strengths 
(Odd > J2«dw)) an d on the range of the interaction, 
which determines the power-law in the generalized ex- 
ponential decay (2/5 > 1/4). The longer range of the 



dipole-dipole interaction makes it more effective for fast 
dephasing in a mesoscopic ensemble [141 ]. In Fig. [2] we 
also observe that the leading-order asymptotics is more 
accurate for D = 1 and 2 than for D — 3. The re- 
sults for D = 3 are significantly improved by including 
the next asymptotic correction. The analytical and nu- 
merical results are then almost indistinguishable on the 
scale shown. Although the cigar-shaped ensemble does 
not have a Gaussian density profile, as assumed in our 
asymptotic analysis, it is possible to qualitatively model 
the data of Fig. [I] for n > 60, if we take the dimension 
of the gas to be nonintegral, in fact D — 1.3. We remark 
that we are showing the behavior of V{t) because we are 
interested in the limit V{t) <C 1. where c/ 2 ' is dominated 
by the two-body interactions. Contributions to the spin 
wave pair correlation function @ from more than two 
excitations can be assessed via formula (0}. 

In conclusion, we have given results for the long-time 
asymptotic dynamics of Rydberg atom spin waves that 
have been prepared via laser excitation of a cold atomic 
ensemble. Our results provide analytical scaling behav- 
iors that are of value in assessing, e.g., the limits in 
speed-up achievable by using Rydberg atom interactions 
for single-photon protocols relevant to quantum informa- 
tion processing [15j . For clouds having a Gaussian den- 
sity profile, the decay follows a generalized exponential in 
time, with a rational power in the exponent that is gov- 
erned by the range of the interaction potential. Numeri- 
cal results for a cigar-shaped cloud are in good agreement 
with a recent experimental measurement of spin wave de- 
phasing in atomic Rb [8j. 
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I. SUPPLEMENTAL MATERIAL: SHORT-TIME BEHAVIOR OF THE ATOM PAIR INTERFERENCE 
FUNCTION FOR DIPOLE-DIPOLE INTERACTION FOR D = 3 

We have calculated the short-time dynamics for D = 3 of the atom pair interference function for the case of isotropic 
dipole-dipole interaction. By employing the method of matched asymptotic expansions (l2T |. we have obtained the 
result 
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(13) 



where 7 = 0.5772 ... is the Euler-Mascheroni constant. 



II. SUPPLEMENTAL MATERIAL: MICROWAVE DRIVEN DIPOLE-DIPOLE INTERACTION 



We have shown in the main text of the paper that, for a mesoscopic atomic ensemble, dipole-dipole interactions 
cause the rapid decay of the spin wave pair correlation function. Achieving the necessary stability to produce this 
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Figure 3. (Color online). Comparison of short-time [Red dashed line from Eq.(13)] and long-time (Green dashed line from the 
main text of the paper) leading asymptotics of |"P3(£)| 2 with the numerical results (Black solid line). The Blue dot-dashed line 
includes the first order long-time asymptotic correction. We consider the target Rydberg level 100s for Rb, and an ensemble 
waist, w — 30/^m. 



interaction with static electic fields is, however, experimentally challenging [l,]. An alternative approach involves 
driving a transition between neighboring Rydberg levels by means of microwave pulses: the mixing of opposite-parity 
Rydberg orbitals provides the atoms with an oscillating dipole moment [l6l |. This technique has recently found 
interesting applications in experiments pj], fl8| - Reference @ elaborates a scheme in which induced dipole-dipole 
interactions in an atomic ensemble may be used to speed up the production of single photons, by orders of magnitude, 
over a non-interacting atomic ensemble. The single photon protocol involves successive Ramsey pulse sequences - each 
of which implies two microwave pulses separated by an interval during which the Rydberg atoms interact. A simple 
analysis of a single Ramsey sequence shows that the amplitude picked up by a double excitation due to interactions 
is given by: \fiv) — > e~ lA ^ t cos(A^ v t)\fii/) . This result differs from the case of dipole-dipole and van der Waals 
interaction, where only a phase shift is introduced by the interaction. The pair correlation function at large time 
may still be calculated from Eq. (5) of the main text, but with the replacement V(t) — > |[1 + V(2t)]. This change 
means that the correlation function relaxes to a finite value smaller than unity (sub-Poissonian statistics) after a 
single pulse; in order to generate true single photons, multiple repetition of the Ramsey sequence are required. In 
Fig. @] we compare the analytical result (lines) for the asymptotic behavior of after a single Ramsey cycle with 
simulations (dots) based on the proposal contained in Ref. ((J. The fast initial decay of g^ is of particular interest, as it 
determines the limiting speed of single photon production. The present theory suggests that this decay is a generalized 
exponential, of the form [r 2 / 5 e -Ro [ B ]( 2T ) 6 ], Although this result is obtained using an isotropic interaction, it is still 
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able to reproduce the main features of the numerical calculation of g^ 2 \ i.e., the steep initial decay to a minimum 
and recovery to a constant asymptotic value. Overall, the qualitative agreement seems good. The discrepancy at 
long times is principally due to the anisotropy of the dipole-dipole interaction, which is included in the numerical 
simulations. 




Figure 4. (Color online) Comparison between (i) the asymptotic behavior of correlations for microwave-induced dipole-dipole 
interactions in a three-dimensional sample having w = 60 /xm and (ii) numerical analysis of the full anisotropic interaction [(J 
(green dots). The correlation function is calculated from the two-body integral using Eq. (5) in the main text. Dashed line: 
leading-order asymptotics; dot-dashed line: including the first asymptotic correction. 
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